November  2002 


3299 


NOTES  AND  CORRESPONDENCE 


Evidence  of  a  Barrier  Layer  in  the  Sulu  and  Celebes  Seas 

Peter  C.  Chu 

Naval  Ocean  Analysis  and  Prediction  Laboratory,  Department  of  Oceanography,  Naval  Postgraduate  School,  Monterey,  California 

Qinyu  Liu  and  Yinglai  Jia 

Physical  Oceanography  Laboratory,  Ocean-Atmosphere  Interaction  and  Climate  Laboratory, 

Ocean  University  of  Qingdao,  Qingdao,  China 

Chenwu  Fan 

Naval  Ocean  Analysis  and  Prediction  Laboratory,  Department  of  Oceanography,  Naval  Postgraduate  School,  Monterey,  California 

13  September  2001  and  4  June  2002 
ABSTRACT 

Variability  in  the  surface  isothermal  and  mixed  layers  of  the  Sulu  and  Celebes  Seas  is  examined  using  the 
conductivity-temperature-depth  data  from  the  Navy’s  Master  Oceanographic  Observational  Data  Set  (MOODS). 

Vertical  gradient  is  calculated  to  determine  isothermal  layer  depth  with  a  criterion  of  0.05°C  m"'  for  temperature 
profiles  and  mixed  layer  depth  with  a  criterion  of  0.015  kg  m  4  for  density  profiles.  When  the  isothermal  layer 
depth  is  larger  than  the  mixed  layer  depth,  the  barrier  layer  occurs.  This  study  shows  that  the  barrier  layer  occurs 
often  in  the  Sulu  and  Celebes  Seas.  In  the  Sulu  Sea,  the  barrier  layer  has  seasonal  variability  with  a  minimum 
occurrence  (38%)  and  a  minimum  thickness  (3  m)  in  May  and  a  maximum  occurrence  (94%)  and  a  maximum 
thickness  (36.5  m)  in  September.  In  the  Celebes  Sea,  the  barrier  layer  thickness  changes  from  a  maximum  (49.7- 
62.0  m)  in  March-April  to  a  minimum  (9.6  m)  in  June.  Possible  mechanisms  responsible  for  the  barrier  layer 
formation  are  discussed.  In  the  Sulu  Sea,  the  barrier  layer  may  be  formed  by  both  rainfall  and  stratification;  in 
the  Celebes  Sea,  a  rain-formed  mechanism  seems  a  major  factor. 


1.  Introduction 

A  vertically  quasi-uniform  layer  of  temperature  (iso¬ 
thermal  layer),  salinity  (isohaline  layer),  and  density 
(mixed  layer)  is  often  observed  in  the  upper  ocean.  Un¬ 
derneath  it,  there  exists  a  layer  with  strong  vertical  gra¬ 
dient  such  as  the  thermocline  (in  temperature),  halocline 
(in  salinity),  and  pycnocline  (in  density).  The  isothermal 
layer  depth  ( HT )  is  not  necessarily  identical  to  the  iso¬ 
haline  layer  depth  ( Hs )  or  the  mixed  layer  depth  (HD). 
There  are  areas  of  the  World  Ocean  where  HT  is  deeper 
than  Hs  or  HD  (Lindstrom  et  al.  1987;  Lukas  and  Lind- 
strom  1991;  Sprintall  and  Tomczak  1992;  You  1995; 
Kara  et  al.  2000).  The  layer  difference  between  halocline 
and  thermocline  is  defined  as  the  barrier  layer  (BL), 
which  has  strong  salinity  stratification  and  weak  (or  neu¬ 
tral)  temperature  stratification  (Fig.  1).  The  barrier  layer 
thickness  (BLT)  is  originally  defined  by  the  difference 


Corresponding  author  address:  Dr.  Peter  C.  Chu,  Department  of 
Oceanography,  Naval  Postgraduate  School,  Monterey,  CA  93943- 
5100. 

E-mail:  chu@nps.navy.mil 


between  HT  and  Hs.  However,  in  recent  studies  BLT  is 
often  referred  to  as  the  difference  between  HT  and  HD 
(e.g.,  Kara  et  al.  2000). 

Less  turbulence  in  the  BL  than  in  the  isohaline  layer 
due  to  strong  salinity  stratification  isolates  the  isohaline 
layer  from  cool  thermocline  water.  This  process  affects 
the  ocean  heat  budget  (Swenson  and  Hansen  1999)  and 
the  heat  exchange  with  the  atmosphere.  Studies  on  BL 
and  related  physical  mechanisms  have  been  devoted  to 
the  open  ocean  such  as  in  the  equatorial  Pacific  with 
BLT  from  20  to  40  m  (e.g.,  Lukas  and  Lindstrom  1991; 
You  1998)  and  the  North  Pacific  with  BLT  larger  than 
50  m  (Kara  et  al.  2000).  However,  little  attention  has 
been  given  to  BL  that  occurs  in  regional  seas.  We  use 
the  conductivity-temperature-depth  (CTD)  data  from 
the  U.S.  Navy’s  Master  Oceanographic  Observational 
Data  Set  (MOODS)  for  the  Sulu  and  Celebes  Seas  to 
demonstrate  the  existence  and  variability  of  the  BL  in 
regional  seas.  The  outline  of  this  paper  is  as  follows.  A 
description  of  the  Sulu/Celebes  Seas  is  given  in  section 
2.  The  MOODS-CTD  data  are  depicted  in  section  3. 
Criteria  for  determining  isothermal  and  mixed  layers  are 
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Fig.  1.  A  diagram  of  isothermal,  mixed,  and  barrier  layers. 


discussed  in  section  4.  A  description  of  BL  in  the  Sulu 
and  Celebes  Seas  is  given  in  sections  5  and  6.  Physical 
processes  for  the  BL  formation  are  depicted  in  section 
7.  In  section  8  we  present  the  conclusions. 

2.  Sulu  and  Celebes  Seas 

The  Sulu  and  Celebes  Seas  are  the  two  major  marginal 
seas  in  the  outer  southeastern  Asia  region  (Fig.  2),  which 
consists  of  an  island  arc  stretching  across  near  5150  km 
along  the  equator  at  about  15°N-11°S  and  94°-141°E, 
and  which  has  two  contrasting  zones  (Arief  1998).  The 
Sulu  Sea  belongs  to  the  western  zone,  where  approxi¬ 
mately  15%-20%  of  the  global  freshwater  discharge 
enters  (Toole  1987)  through  high  annual  rainfall,  be¬ 
tween  2  and  4  m  yr_1  (ASEAN  Subcommittee  on  Cli¬ 
matology  1982).  The  Celebes  Sea  belongs  to  the  eastern 
zone,  which  is  composed  of  deep-basin  chains  of  depth 
1000-4000  m  and  filled  mainly  by  western  Pacific  water 
masses.  This  zone  has  annual  rainfall  less  than  2  m  yr_I 
(ASEAN  Subcommittee  on  Climatology  1982). 

3.  MOODS-CTD  data 

The  temperature  and  salinity  data  for  the  Sulu  Sea 
(8°-12°N,  1 17°-123°E)  and  the  Celebes  Sea  (0°-6°N, 
1 17°-126°E)  acquired  from  MOODS  went  through  qual¬ 
ity  control  procedures  such  as  a  min-max  check  (e.g., 
disregarding  any  temperature  data  less  than  —  2°C  and 
greater  than  35°C),  error  anomaly  check  (e.g.,  rejecting 
temperature  data  deviating  more  than  7°C  from  clima¬ 
tology),  ship-tracking  algorithm  (screening  out  data  with 
obvious  ship  position  errors),  max  number  limit  (lim¬ 
iting  a  maximum  number  of  observations  within  a  spec¬ 
ified  and  rarely  exceeded  space-time  window),  and  bud¬ 
dy  check  (deleting  contradicting  data).  The  climatolog¬ 
ical  data  set  used  for  the  quality  control  is  the  U.S. 
Navy’s  Generalized  Digital  Environmental  Model 
(GDEM)  climatological  temperature  and  salinity  data¬ 
set.  After  quality  control,  there  are  221  (278)  CTD  sta¬ 
tions  in  the  Sulu  (Celebes)  Sea.  The  majority  of  the 


CTDs  were  nominally  capable  of  reaching  a  maximum 
depth  of  2000  m. 

The  temporal  and  spatial  distribution  of  MOODS- 
CTD  data  in  the  Sulu/Celebes  Seas  is  irregular.  Certain 
periods  and  areas  are  very  well  sampled  while  others 
lack  enough  observations  to  gain  any  meaningful  in¬ 
sights.  For  example,  there  is  a  minimum  (maximum) 
number  of  CTD  stations  in  the  Sulu  Sea  in  May  (July). 
There  are  no  observations  in  the  whole  Sulu  (Celebes) 
Sea  in  certain  years  (e.g.,  1942-46,  1952-63  for  the 
Sulu  Sea,  and  1979-88  for  the  Celebes  Sea)  and  there 
are  many  observations  in  other  years  (e.g.,  75  profiles 
in  the  Sulu  Sea  in  1949,  and  52  profiles  in  the  Celebes 
Sea  in  1941).  Spatial  and  temporal  irregularities  along 
with  the  lack  of  data  in  certain  regions  must  be  carefully 
weighted  in  order  to  avoid  statistically  induced  vari¬ 
ability. 

4.  Criteria  for  determining  HT  and  H D 

a.  Existing  criteria 

There  are  two  types  of  criteria,  difference  and  gra¬ 
dient,  for  determining  HT  and  HD  in  the  upper  ocean. 
The  difference  criterion  requires  the  deviation  of  tem¬ 
perature  (density)  from  its  surface  value  to  be  smaller 
than  a  certain  fixed  value.  The  gradient  criterion  requires 
the  vertical  derivative  of  temperature  (density)  to  be 
smaller  than  a  certain  fixed  value. 

1)  Difference  criteria 

The  criterion  for  determining  HT  varies  from  0.5°C 
(Wyrtki  1964;  Monterey  and  Levitus  1997)  to  0.8°C 
(Kara  et  al.  2000).  The  criterion  for  determining  HD  is 
given  by  (Miller  1976;  Spall  1991) 

Act,  =  0.125o-,(0),  (1) 

or  by  (Sprintall  and  Tomczak  1992;  Ohlmann  et  al. 
1996;  Monterey  and  Levitus  1997) 

Act,  =  0.5  °C(da,/dT),  .  (2) 

respectively.  Here  dcrJdT  is  the  thermal  expansion  co¬ 
efficient  evaluated  using  the  surface  values  of  temper¬ 
ature  and  salinity.  The  difference  criterion  (2)  is  based 
on  the  assumptions  that  the  salinity  effect  on  seawater 
expansion  is  negligible  and  that  HD  is  corresponding  to 
the  depth  with  temperature  difference  of  0.5°C  from  the 
surface  (Sprintall  and  Tomczak  1992). 

2)  Gradient  criteria 

Defant  (1961)  was  among  the  first  to  use  the  gradient 
method.  He  uses  a  gradient  of  0.015°C  m_1  to  determine 
Ht  of  the  Atlantic  Ocean.  Bathen  (1972)  chose  0.02°C 
m-1,  and  Lukas  and  Lindstrom  (1991)  used  0.025°C 


da-Jdz  =  0.01  kg  m-4  (4) 
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are  used  here  to  determine  HT  and  HD.  This  is  because 
1)  the  evident  difference  between  the  isothermal  (iso- 
pycnal)  layer  and  thermocline  (pycnocline)  i6  the  ver¬ 
tical  gradient,  and  2)  a  recent  observational  study  has 
shown  that  the  mean  HT  is  less  sensitive  using  gradient 
criteria  than  using  difference  criteria  using  27  station 
data  (You  1995).  The  mean  HT  vary  little  among  cal¬ 
culations  based  on 

dT/dz  =  0.05°C  nr1,  dT/dz  =  0.025°C  nr\ 

dT/dz  =  0.015°C  nr1.  (5) 

c.  Threshold  for  determining  BL 

In  the  CTD  measurements  there  are  generally  two 
thermometers,  whose  accuracy  is  about  0.02°C.  Thus, 
the  error  in  computing  HT  is  around  1  m  using  the  gra¬ 
dient  criterion  (4).  The  density  a,  is  computed  from 
temperature  and  salinity.  The  error  in  computing  HD 
may  be  estimated  by  2  m.  Thus,  the  threshold  for  the 
BL  is  given  as  3  m  in  this  study. 

5.  Existence  of  BL  in  the  Sulu  Sea 

a.  General  statistics 

The  two  depths,  HT  and  HD,  are  obtained  from  221 
CTD  profiles  using  the  gradient  criteria  (4).  Monthly 
CTD  stations  (denoted  by  the  symbol  “o”)  and  the 
stations  with  the  BL  occurrence  {HD  <  HT,  represented 
by  the  symbol  “  +  ”)  in  the  Sulu  Sea  (Fig.  3)  show  a 
rather  frequent  occurrence  of  BL  in  the  Sulu  Sea.  For 
example,  among  16  (6)  CTD  stations  in  September  (No¬ 
vember),  there  are  15  (6)  stations  where  a  BL  occurs. 
The  rate  of  occurrence  reaches  94%  in  September  and 
100%  in  November  (Table  1).  When  BL  occurs,  the  BLT 
is  computed  as  the  difference  between  HT  and  HD.  The 
climatological  characteristics  are  outlined  as  follows:  (i) 
a  BL  occurs  most  frequently  in  September-November 
(fall)  with  a  frequency  of  90%-100%  and  least  fre¬ 
quently  in  April  and  May  (spring)  with  a  frequency  of 
50%;  (ii)  The  mixed  layer  depth  (MLD)  has  a  minimum 
value  of  3  m  in  May  and  a  maximum  value  of  30.5  m 
in  February;  and  (iii)  the  BLT  has  large  values  in  Sep- 
tember-October  (39.7-47.6  m)  and  small  values  of  9- 
11.3  m  in  April  and  May.  It  is  noticed  that  the  sample 
size  is  not  sufficiently  large;  thus  the  statistical  features 
mentioned  here  may  have  bias,  especially  in  May. 

b.  Seasonal  variability 

Four  groups  of  CTD  data  (Fig.  4)  with  relatively  high 
sampling  rates  are  selected  to  represent  winter  (January- 
February  1949,  February-March  1972)  and  summer 
(July-September  1949,  October  1973).  Scatter  diagrams 
of  Ht  and  HD  (Fig.  5)  show  that  the  BL  occurs  quite 
often  in  both  summer  and  winter  and  is  more  evident 
in  the  summer.  For  example,  the  MLD  is  much  less  than 


the  isothermal  depth  (ITD)  at  almost  all  stations  in  Oc¬ 
tober  1973.  The  maximum  thickness  of  the  BL  (80  m) 
is  found  at  Station  6  (see  Fig.  4)  with  a  MLD  of  1 8  m 
and  ITD  of  98  m.  Comparison  among  the  four  periods 
shows  the  BL  has  a  relatively  thin  BLT  in  July-Sep¬ 
tember  1949. 


The  two  depths  ( HT  and  HD)  for  the  Celebes  Sea  are 
obtained  from  processing  179  CTD  profiles  using  the 
gradient  criteria  (5).  Monthly  CTD  stations  (denoted  by 
the  symbol  “o”)  and  the  stations  with  the  BL  occur¬ 
rence  (Hd  <  Ht,  represented  by  the  symbol  “  +  ”)  in 
the  Celebes  Sea  (Fig.  6)  show  a  rather  frequent  occur¬ 
rence  of  BL  in  the  Celebes  Sea.  For  example,  among 
14  CTD  stations  in  December,  there  are  13  stations 
where  BL  occurs.  The  rate  of  occurrence  reaches  93% 
(Table  2).  The  BLT  has  a  maximum  value  of  62.0  m  in 
April  and  a  minimum  value  of  9.6  m  in  June.  The  cli¬ 
matological  features  of  a  BL  are  outlined  as  follows:  (i) 
the  BL  occurs  most  frequently  in  December  with  a  fre¬ 
quency  of  93%  and  least  frequently  in  April  with  a 
frequency  of  36%;  (ii)  the  MLD  has  a  minimum  value 
of  7.3  m  in  March  and  a  maximum  value  of  28.1  m  in 
September;  and  (iii)  the  BLT  has  a  maximum  value  in 
March-April  (49.7-62.0  m)  and  a  minimum  value  of 
9.6  m  in  June.  These  numbers  might  not  have  statistical 
significance  because  of  the  small  sample  size. 


Four  groups  of  CTD  data  (Fig.  7)  with  relatively  high 
sampling  rates  are  selected  to  represent  winter  (January 
1941),  spring  (May-June  1972),  summer  (July  1941), 
and  fall  (October  1978).  Scatter  diagrams  of  HT  and  HD 
(Fig.  8)  show  that  the  BL  occurs  quite  often  in  winter 
and  fall  and  less  often  in  spring.  For  example,  HD  is 
much  less  than  HT  at  almost  all  stations  in  January  1941 . 
At  Station  16,  HD  and  HT  are  25  and  100  m,  respectively, 
with  a  BLT  about  75  m.  Among  the  four  periods,  the 
BL  has  relatively  thin  BLT  in  July  1941  and  May-June 
1972. 

7,  Discussion  of  physical  processes  causing  the  BL 

formation 

a.  Two  mechanisms  for  the  occurrence  of  BL 

Two  major  factors  determining  HT  (Hp)  are  1)  surface 
winds  and  net  heat  (buoyancy)  flux  and  2)  thermal  (den¬ 
sity)  stratification  underneath  the  isothermal  (mixed) 
layer.  Thus,  the  condition  (HT  >  HD)  is  caused  by  the 
surface  freshwater  flux  (precipitation  excess  over  evap¬ 
oration)  and  strong  salinity  stratification  underneath  the 
mixed  layer.  Thus,  the  BL  can  be  classified  as  (i)  rain 


6.  Existence  of  BL  in  the  Celebes  Sea 

a.  General  statistics 


b.  Seasonal  variability 


Latitude  Latitude  Latitude  Latitude 
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Fig.  4.  CTD  stations  in  the  Sulu  Sea  for  the  four  selected  periods:  (a)  Jan-Feb 
1949,  (b)  Feb-Mar  1972,  (c)  Jul-Sep  1949,  and  (d)  Oct-Dee  1973. 

formed,  and  (ii)  stratification  formed  (weaker  thermal 
and  strong  salinity  stratification). 

1)  Rain-formed  mechanism 

There  are  two  regimes  in  the  ocean  mixed  layer  dy¬ 
namics  (e.g.,  Chu  et  al.  1990;  Chu  and  Garwood  1991; 

Chu  1993):  entrainment  and  detrainment.  The  detrain- 
ment  process  occurs  with  weak  winds  and  strong  surface 
warming  (or  excess  precipitation  over  evaporation).  The 
two  depths  Ht  and  HD  are  determined  by 


Table  1.  Monthly  occurrence  frequency  and  climate  features  of  the  barrier  layer  in  the  Sulu  Sea. 


Month 

Jan 

Feb 

Mar 

Apr 

May 

Jun 

Jul 

Aug 

Sep 

Oct 

Nov 

Dec 

Total 

N 

12 

13 

19 

8 

2 

61 

16 

19 

16 

31 

6 

18 

221 

B 

8 

11 

17 

4 

1 

48 

11 

16 

15 

28 

6 

14 

179 

BIN  (%) 

67 

85 

89 

50 

50 

79 

69 

84 

94 

90 

100 

78 

Mean  MLD  (m) 

29 

30.5 

19.5 

16.5 

3 

16.4 

14.2 

12.9 

11.1 

14.9 

21.2 

12.1 

Mean  BLT  (m) 

27.5 

29.9 

14.1 

11.3 

9.0 

12.2 

28.7 

29.3 

39.7 

47.6 

33.7 

45.3 

'  €£*  0  ; 

C.ul  C,ul 

Ht  =  — - —J - ,  Hd  =  (6) 

C2agF/p0cp  C2B 

with  C,  and  C2  being  tuning  coefficients,  g  the  gravity, 
a  the  thermal  expansion  coefficient,  p0  the  characteristic 
density,  cp  the  specific  heat  under  constant  pressure,  and 
F  the  net  surface  heat  flux  (downward  positive).  In  a 
detrainment  regime,  F  >  0.  The  ratio  between  the  two 
depths  is  calculated  by 

H1=l  +  P(P-E)S 
Hd  aF/Pocp  ’ 
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Fig.  5.  Scatter  diagrams  (HT  vs  Hu)  of  the  Sulu  Sea  for  the  four  selected  periods:  (a)  Jan-Feb  1949,  (b)  Feb-Mar 

1972,  (c)  Jul-Sep  1949,  and  (d)  Oct-Dec  1973. 


where  P  and  E  are  the  precipitation  and  surface  evap¬ 
oration  rate  and  /3  is  the  salinity  contraction  coefficient. 
Thus,  the  rain-formed  mechanism  (P  —  E>  0)  becomes 
evident  in  the  detrainment  regime  with  weak  winds, 
weak  surface  warming  (low  positive  value  of  F),  and 
strong  surface  freshwater  flux  [large  positive  value  of 

C p  -  E)]. 

2)  Stratification-formed  mechanism 

The  stratification-formed  mechanism  is  evident  in  the 
entrainment  regime.  Suppose  that  the  rain-formed  mech¬ 
anism  is  absent  (P  =  E)  and  that  the  initial  HT  coincides 
with  Hd: 

tfrU  =  Hd\i=0.  (8) 

The  ratio  of  the  initial  entrainment  velocities  ( we  =  dEH 
dt )  between  isothermal  and  mixed  layers  is  given  by 
(Chu  1993): 


=  [~ (P  ~  P-HB)1  (9) 

wp  a(T  -  T_Ht)  • 

A  larger  density  jump 

[~(P  -  P-hd)]  >  a(T  -  T_Ht)  (10) 

leads  to  a  larger  entrainment  velocity  for  the  isothermal 
layer  than  for  the  mixed  layer.  Since  the  density  strat¬ 
ification  is  determined  by  the  temperature  and  salinity 
stratifications,  the  stratification-formed  mechanism  be¬ 
comes  important  when  the  salinity  stratification  is 
strong.  Thus,  the  stratification-formed  mechanism  be¬ 
comes  evident  in  the  entrainment  regime  with  strong 
winds,  strong  surface  cooling  (negative  value  of  F ),  and 
strong  salinity  stratification.  Usually,  both  rain-formed 
and  stratification-formed  mechanisms  may  take  place  at 
the  same  time,  and  precipitation  may  also  strengthen  the 
stratification. 


3306 


JOURNAL  OF  PHYSICAL  OCEANOGRAPHY 


Volume  32 


1 1 8  E  1 20  E  122°E  124°E 

Longitude 


118  E  120  E  122  E  124°E 

Longitude 


118  E  120  E  122  E  124°E 

Longitude 


Fig.  6.  Monthly  CTD  stations  (denoted  by  “o”)  and  the  stations  with  BL  occurrence  (HD  <  H,\  denoted  by  “  +  ”) 

in  the  Celebes  Sea. 


b.  Sulu  Sea 

The  Comprehensive  Ocean-Atmosphere  Data  Set 
(COADS)  provides  monthly  mean  sea  surface  wind 
stress  and  net  heat  and  freshwater  fluxes  (Fig.  9).  The 
surface  wind  stress,  net  heat  flux  ( F ),  and  net  freshwater 
flux  ( P  —  E)  have  seasonal  variability  in  the  Sulu  Sea. 
The  wind  stress  strengthens  in  December-February 
[6.57-6.67  (X10-2  N  m~2)]  and  weakens  in  April-No- 
vember  [0.63-3.66  (X10-2  N  m-2)].  The  net  heat  flux 
is  negative  in  November-January  (  —  25.4  to  0  W  m-2, 
surface  cooling)  and  becomes  large  positive  in  March- 
May  (60-110.3  W  m-2,  strong  surface  warming).  The 
surface  freshwater  flux  ( P  —  E)  is  weak  in  January- 
May  (<0.6  m  yr-1)  and  strong  in  June-November  (>1 
m  yr-1)- 

In  April  and  May  the  Sulu  Sea  is  under  the  influence 


of  weak  winds  [1.3-3  (X  10~2  N  m-2)],  strong  net  heat 
flux  (60-110.3  W  m-2),  and  weak  surface  freshwater 
flux  (0.52-0.56  m  yr-1).  The  strong  heat  flux  and  weak 
winds  cause  a  shallow  isothermal  layer  (small  HT).  The 
weak  freshwater  flux  cannot  strengthen  the  salinity  strat¬ 
ification.  Thus,  the  BL  occurs  least  frequently  in  April 
and  May. 

In  September  and  October  the  Sulu  Sea  is  under  in¬ 
fluence  of  weak  winds  [0.63-2.19  (X10-2  N  m-2)], 
weak  net  heat  flux  (24.6-29.2  W  m-2),  and  strong  sur¬ 
face  freshwater  flux  (3.1  m  yr-1  in  September).  These 
conditions  lead  to  the  detrainment  regime  for  the  upper 
ocean  and  favor  the  rain-formed  mechanism  (excess  pre¬ 
cipitation  over  evaporation). 

In  December  the  Sulu  Sea  is  under  influence  of  strong 
winds  (6.57  X  10-2  N  m-2),  surface  cooling  (  —  25.4  W 
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Fio.  7.  CTD  stations  in  the  Celebes  Sea  for  the  four  selected  periods:  (a)  Jan-Feb 
1949,  (b)  Feb-Mar  1972,  (c)  Jul-Sep  1949,  and  (d)  Oct-Dee  1973. 


m-2),  and  relatively  weak  surface  freshwater  flux  (0.8 
m  yr-1).  These  conditions  favor  the  entrainment  regime 
in  the  upper  ocean  and  the  stratification-formed  mech¬ 
anism  for  the  occurrence  of  a  BL. 

c.  Celebes  Sea 

The  surface  wind  stress,  net  heat  flux  (F),  and  net 
freshwater  flux  ( P  —  E)  have  less  seasonal  variability 
in  the  Celebes  Sea  than  in  the  Sulu  Sea  (Fig.  9).  The 
wind  stress  varies  between  0.45  X  10~2  (Nov)  and  2.6 
X  10-2  N  m  2  (Feb).  The  surface  net  heat  and  fresh¬ 
water  fluxes  are  positive  year-round  with  F  varying  from 
17.5  (Jan)  to  87.4  W  m-2  (Mar)  and  ( P  —  E)  varying 
from  0.17  (Apr)  to  2.1  m  yr-1  (Jan).  Thus,  the  Celebes 
Sea  is  under  the  influence  of  weak  winds,  surface  warm¬ 
ing,  and  surface  freshwater  flux  all  year.  These  condi¬ 
tions  favor  the  rain-formed  mechanism. 


In  December-January  the  Celebes  Sea  experiences 
weak  net  heat  flux  (17.5-24.5  W  m-2)  and  strong  sur¬ 
face  freshwater  flux  (1.14-2.06  m  yr-1).  These  condi¬ 
tions  may  cause  the  detrainment  regime  for  the  upper 
ocean  and  favor  the  rain-formed  mechanism  (excess  pre¬ 
cipitation  over  evaporation). 

In  April,  the  Celebes  Sea  experiences  strong  net  heat 
flux  (80.4  W  m-2)  and  weak  surface  freshwater  flux 
(0.17  m  yr-1).  The  strong  heat  flux  and  weak  winds 
cause  a  shallow  isothermal  layer  (small  HT).  The  weak 
freshwater  flux  strengthens  the  salinity  stratification 
slightly.  Thus,  the  BL  occurs  least  frequently  in  April. 

8.  Summary 

1)  Atmospheric  forcing  (rain-formed  mechanism)  and 
ocean  physical  processes  (stratification-formed 


Table  2.  Monthly  occurrence  frequency  and  climate  features  of  the  barrier  layer  in  the  Celebes  Sea. 


Month 

Jan 

Feb 

Mar 

Apr 

May 

Jun 

Jul 

Aug 

Sep 

Oct 

Nov 

Dec 

Total 

N 

22 

10 

10 

22 

0 

13 

77 

50 

9 

21 

30 

14 

278 

B 

19 

9 

7 

8 

0 

9 

68 

41 

7 

17 

26 

13 

224 

BIN  (%) 

86 

90 

70 

36 

69 

88 

82 

78 

81 

87 

93 

Mean  MLD  (m) 

21.2 

12.3 

7.3 

13.0 

14.4 

15.3 

24.3 

28.1 

16.1 

15.8 

16.2 

Mean  BLT  (m) 

41.1 

26 

49.7 

62.0 

9.6 

27.9 

30 

11.3 

20.7 

15.3 

16.8 
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Fig.  8.  Scatter  diagrams  ( HT  vs  HD)  of  the  Celebes  Sea  for  the  four  selected  periods:  (a)  Jan-Feb  1949,  (b)  Feb- 
Mar  1972,  (c)  Jul-Sep  1949,  and  (d)  Oct-Dec  1973. 


mechanism)  may  cause  the  occurrence  of  a  barrier  imum  value  of  7.3  m  in  March  and  a  maximum  value 

layer  (deeper  thermocline  than  pycnocline)  in  the  of  28.1  m  in  September.  The  barrier  layer  thickness 

Sulu/Celebes  Seas.  Using  the  U.S.  Navy’s  Master  has  large  values  in  March-April  (49.7-62.0  m)  and 

Oceanographic  Observational  Data  Set,  we  have  small  value  of  9.6  m  in  June.  Relatively  higher  pre¬ 
shown  barrier  layer  thickness  up  to  47.6  m  (62.0  m)  cipitation  and  lower  wind  stress  than  for  the  Sulu 

in  the  Sulu  (Celebes)  Sea.  Sea  are  found  in  winter  (see  Fig.  9). 

2)  The  barrier  layer  appears  quite  often  in  both  Sulu  5)  In  the  Sulu  Sea,  both  rain-formed  and  stratification- 

and  Celebes  Seas.  The  estimated  occurrence  rate  formed  mechanisms  are  responsible  for  the  barrier 

varies  from  50%  (April  and  May)  to  100%  (Novem-  layer  formation.  However,  in  the  Celebes  Sea  the 

ber)  in  the  Sulu  Sea  and  from  36%  (April)  to  93%  rain-formed  mechanism  is  responsible  for  the  barrier 

(December)  in  the  Celebes  Sea.  layer  formation. 

3)  In  the  Sulu  Sea,  the  mixed  layer  depth  has  a  mini¬ 
mum  value  of  3  m  in  May  and  a  maximum  value  of 

30.5  m  in  February.  The  barrier  layer  thickness  has  Acknowledgments.  This  work  was  funded  by  the  Na- 
large  values  in  September-November  (39.7-47.6  m)  val  Oceanographic  Office,  the  Office  of  Naval  Research, 
and  small  values  of  9-1 1.3  m  in  April  and  May.  and  the  Naval  Postgraduate  School.  Liu  was  also  funded 

4)  In  the  Celebes  Sea,  the  mixed  layer  depth  has  a  min-  by  the  Chinese  National  Key  Program  for  Developing 
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Fig.  9.  Monthly  mean  CO  ADS  (a)  surface  wind  stress  (N  m-2), 
(b)  P  —  E  (m  yr~'),  and  (c)  net  heat  flux  (W  m-2)  averaged  over 
the  Sulu  Sea  (solid  curve)  and  the  Celebes  Sea  (dashed  curve). 


Basic  Science  (G1999043807)  and  Chinese  Special  Pro¬ 
gram  “2001DIA50041”  of  the  National  Ministry  of  Sci¬ 
ence  and  Technology. 
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